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Abstract: We use gauge-gravity duality to study the stability of zero-temperature, finite 
baryon density states of = 4 supersymmetric SU (Nc) Yang-Mills theory coupled to a 
single massive fundamental-representation Af = 2 hypermultiplet in the large- A^c and large- 
coupling limits. In particular, we study the spectrum of mesons. The dual description is 
a probe D7-brane in anti-de Sitter space with a particular configuration of worldvolume 
fields. The meson spectrum is dual to the spectrum of fluctuations of worldvolume fields 
about that configuration. We use a combination of analytical and numerical techniques 
to compute the spectrum, including a special numerical technique designed to deal with 
singular points in the fluctuations' equations of motion. Despite circumstantial evidence 
that the system might be unstable, such as a finite entropy density at zero temperature 
and the existence of instabilities in similar theories, we find no evidence of any instabilities, 
at least for the ranges of frequency and momenta that we consider. We discover a pole on 
the imaginary frequency axis in a scalar meson two-point function, similar to the diffusive 
mode in the two-point function of a conserved charge. 
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1. Introduction 

Many real physical systems involve a finite chemical potential and strongly-interacting 
degrees of freedom. Examples include Quantum Chromodynamics (QCD) with a finite 
baryon number chemical potential and many condensed matter systems. 

Such systems often exhibit unusual transport properties. For example, the quark-gluon 
plasma (QGP) created at the Relativistic Heavy-Ion Collider (RHIC) appears to have the 
lowest viscosity of any known substance [1]. Many condensed matter systems, such as 
the heavy-fermion compounds and the high- Tc cuprates, have an electrical resistivity that 
scales linearly with temperature, in stark contrast to a typical Fermi liquid. Such behavior 
appears to be intimately related to quantum criticality, i.e. the presence of a continuous 
phase transition at zero temperature as a function of some non-thermal control parameter 
(pressure, magnetic field, etc.) [2-4]. The low-energy dynamics at a quantum critical point 
is generically described by some strongly-coupled, scale-invariant theory. That theory can 
also describe finite-temperature physics, as long as temperature is the largest scale, and 
hence can influence the finite-temperature resistivity. 

Few reliable methods exist to compute physical observables for these systems, with 
real-time physics being particularly difficult to study. Lattice simulation is not reliable 
because either a finite chemical potential or Lorentzian signature, as required for real-time 
transport, generically introduces the "sign problem," in which a complex phase renders 
Monte-Carlo sampling practically impossible. 
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An alternative method is the anti-de Sitter/Conformal Field Theory (AdS/CFT) cor- 
respondence, or more generally gauge-gravity duality [5-7]. Gauge-gravity duality is a holo- 
graphic duality between a strongly-coupled non-Abelian gauge theory and a weakly-coupled 
theory of gravity in one higher dimension. The field theory "lives" on the boundary of the 
bulk holographic spacetime, for example if the bulk spacetime is AdS then the dual CFT 
lives on the AdS boundary. The field theory energy scale is dual to the radial/holographic 
direction. Every gauge-invariant field theory operator is dual to a field in the bulk. For 
example, the stress-energy tensor is dual to the metric, while the conserved current of a 
global U{1) symmetry is dual to a bulk U{1) gauge field. A thermal equilibrium state 
with finite entropy appears in the gravity theory as a black hole spacetime [8] . In the field 
theory, perturbations away from thermal equilibrium are described by a low-cncrgy effec- 
tive theory called hydrodynamics. Via gauge-gravity duality, the study of hydrodynamics 
becomes the study of small perturbations of black holes [9]. 

Gauge-gravity duality currently does not describe any real system, but is valuable nev- 
ertheless as one of the few ways to answer questions about finite-density, real-time physics 
in strongly-coupled systems. Moreover, gauge-gravity duality has revealed remarkable uni- 
versality in the transport properties of such systems: any theory with a gravitational dual 
(with Einstein-Hilbert action, in a state with S0{2) rotational symmetry [10]) has the 
same, very small, ratio of shear viscosity r] to entropy density s, namely r]/s = 1/4tt [11]. 
That value is remarkably close to the value estimated for the QGP created at RHIC [12]. 

Holographic models fall into two categories, "top-down" and "bottom-up." The former 
begin with a genuine string theory or supergravity action. In such cases the dual field theory 
can (usually) be determined precisely, allowing not only for specific translation between field 
theory and gravity, but also for comparison to weak-coupling results. The price is a large 
number of fields, often with complicated couplings. Bottom-up constructions begin with 
a simple bulk action, including few fields with simple couplings, to capture the essential 
physics. That provides technical simplicity without sacrificing universality. Whether a 
dual field theory exists or is well-defined is not guaranteed, however, so any results must 
be handled with care. 

The simplest bottom-up holographic model of strongly-interacting matter is Einstein- 
Maxwell theory (with negative cosmological constant) [13]. The dual is some CFT with a 
global U{1) symmetry. The AdS-Reissner-Nordstrom (AdS-RN) charged black hole solu- 
tion describes thermal equilibrium states with a finite U{1) charge density. The extremal 
limit of AdS-RN is dual to the field theory's zero-temperature limit at finite density, and 
is the paradigm of holographic quantum critical matter [14, 15]. 

Extremal AdS-RN has a number of unusual properties. Chief among them is the 
presence of an extremal horizon, indicating that the entropy density does not vanish in 
the zero-temperature limit. In other words, the zero-temperature, finite-density state has 
a large degeneracy of microstates. That renders the system vulnerable to instabilities: 
generically, any small perturbation will break the degeneracy, and the system will settle 
into a new, presumably non-degenerate, ground state. 

On the gravity side, the instability is intimately related to the presence of a (1 -|- 1)- 
dimensional AdS factor, AdS2, in the near-horizon region [16]. A condition for the stability 
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of a field in AdS is the Breitenlohner-Frccdman (BF) bound [17], a bound on how negative 
the mass (squared) of a bulk field can be. A bulk field may satisfy the BF bound in the 
near-boundary region but violate the AdS2 BF bound. The exact nature of the instability 
depends on the bulk field content. For example, a charged scalar field may become unstable 
and develop a non-trivial profile: the horizon is replaced by "scalar hair" [18]. The dual 
field theory state is then a zero-temperature superfluid in which the global U{1) symmetry 
is spontaneously broken. The near-horizon AdS2 factor indicates that in the field theory 
some (0-f l)-dimensional CFT emerges at low energies. 

Extremal AdS-RN also has unusual transport properties. Holographic calculations 
revealed that the dual theory has a finite shear viscosity, and that indeed rj/s = I/Att [19]. 
The dual theory also admits a sound mode with speed v identical to the value dictated by 
scale invariance, = 1/d in d spatial dimensions [20]. 

These results are mysterious. On general grounds we do not expect hydrodynamics to 
remain a reliable effective description at arbitrarily low temperatures. Hydrodynamics is 
an effective description in which we expand expectation values of conserved quantities, like 
the stress-energy tensor, in a small parameter, momentum times the mean free path. In the 
CFT the only meaningful parameter is the ratio of temperature to chemical potential, T/^u. 
In the limit T/n ^ 1, we expect the mean free path to scale as 1/T plus corrections in /u/T. 
As we decrease T/fi, however, the mean free path grows, so we expect the hydrodynamic 
expansion to break down at some point. If we continue to use hydrodynamics beyond 
that point, we may see sick behavior. For example, if r] remains finite as T/|Lt — > then 
the divergence of the entropy current grows without bound [21]. The holographic results 
for zero-temperature transport suggest that perhaps some effective theory, very much like 
hydrodynamics, exists in the limit T/fi<^ 1. What theory replaces hydrodynamics at low 
temperatures? For an attempt to answer this question via a holographic version of the 
Wilsonian renormalization procedure, sec rcf. [21]. 

Our goal is to perform a complete stability analysis in a top-down model of holographic 
matter at zero temperature. AdS-RN may be extracted from top-down holographic models 
(see for example refs. [22-25] and references therein), which usually involve many fields 
with complicated couplings. To our knowledge no complete stability analysis, including all 
fields, has been performed in any top-down model. ^ We will therefore employ a tractable 
substitute: a probe D-brane. 

We will study type IIB supergravity in the near-horizon geometry of very many D3- 
branes, AdS^ x S^, with a single probe D7-brane extended along AdS^ x [26]. The 
dual field theory is M = 4 supersymmetric SU{Nc) Yang-Mills (SYM) in the large-iVc 
and large-coupling limits, coupled to a single J\f = 2 hypermultiplet in the fundamental 
representation of SU{Nc), i.e. flavor fields. The field theory has a global U{1) flavor 
symmetry analogous to baryon number in QCD, which we will denote U{1)b. We will 
introduce a baryon number chemical potential and an M = 2 supersymmetric mass for 

^Many top-down models rely on (consistent) truncations, in which the full set of supergravity fields is 
restricted to some subset of fields that only source fields within that subset. Obviously, only instabilities 
of fields within that subset can then be detected, that is, the truncation leaves us blind to instabilities of 
other fields, and so we cannot determine the true ground state with certainty. 
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the hypermultiplet. The baryon number current and the flavor mass operator are dual to 
the U{1) gauge field and a scalar field on the D7-brane worldvolumc [27]. The D7-brane 
action, including the Dirac-Born-Infeld (DBI) and Wess-Zumino (WZ) terms, describes 
the dynamics of these fields. Fortunately, exact solutions for these fields are known that 
describe field theory states with a finite baryon density at zero temperature, when the 
chemical potential is larger than the mass [28]. We will study the spectrum of linearized 
fluctuations of worldvolume fields about these background solutions, dual to the spectrum 
of mesonic operators [29]. 

Our system shares many features with extremal AdS-RN, in particular a finite entropy 
density and a sound mode with speed identical to the conformal value [30-32]. An AdS2 
region is also present in some sense, since fluctuations of worldvolume fields effectively "see" 
an AdS2 region deep in the bulk, due to their couplings to the background worldvolume 
fields [21,33-35], as we will review in detail below. ^ 

We have four main reasons to expect instabilities in our system. The first is the finite 
entropy density at zero temperature. The second is the presence of charged scalars in the 
field theory susceptible to Bose-Einstein condensation. The third is the presence of charged 
fermions: our intuition from large- AT^ QCD with an asymptotically large chemical potential 
is that quarks may exhibit an instability towards an inhomogeneous ground state, the so- 
called "chiral density wave" [36, 37]. The fourth reason comes from other holographic 
systems, where either bulk Chern-Simons or axion-like couplings triggered instabilities 
towards inhomogeneous ground states [38-43]. The WZ terms in the D7-brane action 
indeed act as Chern-Simons couplings for the worldvolume gauge field. 

Our analysis is complete in that we include fluctuations of all bosonic D7-brane fields.^ 
The only other fields in our system are those of supergravity, but in the probe limit these 
are insensitive to the D7-brane. We are thus including all relevant bosonic fields in our top- 
down model. Crucially, recall that none of the D7-brane worldvolume fields are charged 
under the worldvolume U{1), or in field theory language mesons are not charged under 
baryon number, so we know a priori that we will not see a superfluid instability. We 
include in our analysis fluctuations with flnite momentum, in anticipation of instabilities 
towards inhomogeneous ground states. Our analysis is incomplete in that we do not study 
arbitrarily high momenta, although that is a limitation of our numerical procedure rather 
than a matter of principle. 

Despite all the circumstantial evidence, we find no signs of any instabilities. This top- 
down model appears to be perturbatively stable, at least over a wide range of momenta. 
We hasten to add that this low-momentum stability is probably special to the large- Ac and 
probe limits, as we discuss in what follows. 

In our analysis we recover the "zero sound" mode of refs. [30,31], appearing as a pole 

^If we leave the probe limit, and allow the D7-brane to back-react on the fields of supergravity, then we 

might expect to find an extremal black hole, which on general grounds we expect to have a finite entropy 
and near-horizon AdS2 region. The features we see for the probe D7-brane may simply be the probe-limit 
remnants of those for the extremal black hole. 

^Fermionic fields also live on the D7-brane worldvolume [44,45], but on general grounds we do not expect 
them to become unstable, as explained for example in ref [16], so we omit them from our analysis. 
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in the U{1)b charge density two-point function. We also discover a new mode with purely 
imaginary dispersion in the two-point function of a certain scalar operator that is a vector of 
the SU{2)]i R-symmetry. We cannot resist drawing an analogy between this mode and the 
"spin diffusion" mode in some itinerant electron systems. Here we adopt the perspective of, 
for example, refs. [46,47]. In a system of electrons, at low energies spin-orbit interactions are 
suppressed and spin rotations decouple from spacetime rotations. Spin essentially becomes 
an internal SU{2) quantum number, that is, a global symmetry, for which we may ask 
about the rate of diffusion. In our case, the SU (2) is part of the R-symmetry, so we dub 
the mode we find "R-spin diffusion." The presence of the zero sound and R-spin diffusion 
lend further evidence for the existence of an effective, hydrodynamic-like theory describing 
holographic quantum critical matter. 

The absence of instabilities may be consistent with the interpretation, proposed in 
refs. [28,48], that the solution of rcf. [28] (about which we perturb) represents scalars 
that have already Bose-Einstein condensed. The argument proceeds as follows. Recall 
that an A/" = 2 hypermultiplet includes a Dirac fermion and two complex scalars, all 
with the same U{1)b charge. In our theory the scalars have quartic self- interactions as 
well as Yukawa and other couplings.^ We then introduce a U{1)b chemical potential, 
which we begin to increase. When the chemical potential is less than the flavor mass, the 
ground state is indistinguishable from the vacuum: all operators have zero expectation 
value. Once the chemical potential equals the flavor mass, the scalars should condense, 
producing a nonzero density. Indeed, holographic calculations reveal a phase transition of 
precisely that character [28]. In addition, the M = 2 supersymmetric completion of the 
hypermultiplet Dirac mass operator acquires an expectation value identical in form to that 
of condensed scalars with quartic self-interaction. The fact that we find no instabilities 
may be because the instability already occurred: the scalars condensed, and the system is 
already in the ground state. Such an interpretation has several problems, however, as we 
review in detail below. Chief among them is the fact that the scalars are charged under 
SU{Nc) X SU{2)r X U{1)b, so if they condense we expect these symmetries to be broken 
to some subgroup. That does not occur for the solution of ref. [28]. Our calculation serves 
to characterize further whatever state the solution of rcf. [28] represents: the low-energy 
spectrum includes not only the zero sound mode but also the R-spin diffusion mode, and 
at very low energies a (0-1-1 )-dimensional CFT emerges. 

The time component of the extremal AdS-RN metric has a double zero. That gener- 
ically introduces an irregular singular point in the equations of motion of fluctuations, 
making numerical solutions difficult to obtain. Two techniques designed to address this 
problem are matched asymptotic expansions, which is reliable for sufficiently low frequen- 
cies [16], and Leaver's method [50,51]. In our system the fluctuation equations also involve 
irregular singular points. We employ matched asymptotic expansions, but Leaver's method 
is difficult to use, as we explain below. We thus develop a new numerical technique to com- 
pute the spectrum, which we call the "zig-zag" method (for reasons that will become clear), 
which may be applicable more generally. 

*The Lagrangian of our theory is written explicitly in ref. [49]. 
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This paper is organized as follows. In section ^ we briefly review the field theory and 
its holographic dual, the circumstantial evidence for instabilities, and the interpretation of 
the solution of ref. [28] as a condensate of scalars. In section ^ we derive the equations of 
motion for fluctuations of D7-brane fields. In section ^ we perform matched asymptotic 
expansions, explain the zig-zag method, and present our numerical results for the spectrum. 
We conclude with some discussion in section |5|. 

2. The Theory and Its Holographic Dual 

We study a strongly-coupled CFT, namely = 4 SYM theory with gauge group SU{Nc) 
in the large- A^c limit and with large 't Hooft coupling A = gyj^Nc ^ 1, where qym is 
the SYM coupling. We introduce a single M = 2 supersymmetric hypermultiplet in the 
fundamental representation of the gauge group, which in analogy with QCD we will call a 
flavor field. The field content of an A/" = 2 hypermultiplet is a Dirac fermion ip and two 
complex scalars, q and q. We will call the fermion a quark and the scalars squarks. 

We work in the probe limit: with Nf flavors of hypermultiplets, the probe limit consists 
of keeping Nf fixed as A'^c — ^ cO; so that Nf/Nc <C 1. If we work to leading order in that 
small parameter, then we may neglect quantum effects due to the fiavor fields, such as the 
running of the coupling. For example, at weak coupling, A <C 1, the probe limit consists of 
discarding all diagrams with fiavor fields in any loops. 

AA = 4 SYM theory has an SO{6)r R-symmetry. The J\f = 2 supersymmetric flavor 
fields have superpotential couplings that break the symmetry to 50(4) xU{1)r ~ SU (2) l x 
SU{2)r X U{1)r, of which the SU{2)r x U{1)r acts as the AA = 2 R-symmetry. The U{1)r 
acts as a chiral symmetry for the quarks, while the squarks are neutral. 

We will introduce an AA = 2 supersymmetry-preserving mass M for the flavor flelds, 
which explicitly breaks the chiral U{1)r, much like QCD. We will denote the J\f = 2 
supersymmetric completion of the Dirac fermion mass operator as Om- Of the three 
complex adjoint scalars of A/" = 4 SYM, only one is charged under U{1)r. Denoting that 
scalar as ^, Om schematically takes the form (the precise operator is written in ref. [27]) 

Om = ii^i^ + q^M + ^')g + q^{M + ^)q. (2.1) 

The theory has a C/(1)b fiavor symmetry, analogous to baryon number in QCD. We will 
introduce a baryon number chemical potential /i. Although we will mostly be interested 
in the theory at zero temperature, we will sometimes discuss the theory in a thermal 
equilibrium state with temperature T. 

To obtain the holographic dual of the above theory we begin in type IIB string theory 
with the following supersymmetric intersection of Nc D3-branes with a single D7-brane, 
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Open strings with both ends on the D3-branes give rise at low energies to AA = 4 SYM, 
while open strings with one end on the D3-branes and one on the D7-branes give rise to the 
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single N = 2 hypermultiplet. If we separate the D3- and D7-branes by a distance L in the 
(ccgjXg) plane, then the open strings between them acquire a length and hence the flavor 
fields acquire a mass that is L times the string tension, M = L/ (27ra'). The relative angle 
in the (xgjXg) plane represents the phase of the mass. The U(1)r symmetry corresponds 
simply to rotations in the (xg, 0:9) plane. 

Taking the usual Maldacena limits, we obtain the holographic dual of the large-A^c) 
large-coupling J\f = A SYM theory, type IIB supcrgravity in the near-horizon geometry of 
the D3-branes, AdS^ x with N^. units of Ramond-Ramond (RR) five-form flux on the 
S^. We will write the metric and five-form as 

ds"^ = Gavdx^dx'' = ^{-dt^ + dx^) + ^idr"^ + r'^dnh, 

= -^Ar^dx^ A dx^ A dx'^ A dx^ A dr) - AR^n^, (2.3) 
R 

where R is the 74^5*5 radius of curvature, r is the AdS^ radial coordinate, with the boundary 
at r — >^ 00, the field theory Minkowski coordinates are t and x, and dil^ is the volume form 
of S^. We will use an metric of the form 

dnl = de^ + sin^ ed(f)^ + cos^ ednl (2.4) 

where dQ^ is the volume form of S^. The RR four-form potential C^^^ giving rise to the 
five-form via F^^^ = dC^^^ then takes a simple form: 

C(^) = -^(r^da;° A dx^ A dx^ A dx^ + i?^ cos^ 9d(j) A dfls). (2.5) 
R 

The radius of curvature is fixed to be R^ = AirggNca'^ where gs is the string coupling. 
In the AdS/CFT dictionary Airgs = gyu^ ^° large 't Hooft-coupling limit of the CFT 
corresponds to R'^ /a'"^ ^ 1, the limit in which excited string states become heavy relative 
to the AdS^ scale. ^ The AdS/CFT dictionary also equates the isomctries of the bulk 
spacetime with the symmetries of the CFT: the isometry group of AdS^ is dual to the 
(3-|-l)-dimensional conformal group, and the 50(6) isometry of the is dual to the R- 
symmetry of A/" = 4 SYM. 

In the near-horizon limit, the D7-brane is extended along AdS^ x S'^ inside AdS^ x . 
The action for a D7-brane in this background includes the DBI and WZ terms, 

Sd7 = -Td7 J d^^V-det(5ab + {27Ta')Fai>) + ^^^Tb7 J A F A F, (2.6) 

where Td7 = ^^{2^)T is the D7-brane tension, (^^, . . . , ^^) are the worldvolume coordinates, 
gab is the induced metric of the D7-brane, Fab = da^b ~ db-^a is the field strength associated 
with the worldvolume f7(l) gauge field Aa-, and P[C(^-'] is the puUback of the RR four-form 
to the D7-brane worldvolume. 

The statement that the D7-brane is dual to the N = 2 hypermultiplet means that every 
field on the D7-brane is dual to some mesonic operator, that is, a gauge-invariant operator 

^Starting now, we will use units in which R = \. In these units, a' = \~^^^ is dimensionless. 
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bilinear in hypermultiplet fields. The fields living on the D7-brane are two scalars, repre- 
senting fluctuations of the D7-brane in the two transverse directions, and the eight compo- 
nents of the U{1) gauge field. These are dual to spin zero and spin one mesons. Fermionic 
fields also live on the D7-brane worldvolume, which are dual to fermionic mesons [44,45], 
but we will leave an analysis of their fluctuations for future work. Generically, higher-spin 
mesons are dual to heavy string excitations accessible only outside of the supergravity 
limit. 

Consider for example Om, which is dual to a D7-brane worldvolume scalar. Indeed, 
before the near-horizon limit that scalar represents the separation of the D3- and D7- 
branes in the overall transverse plane. To describe that here, we rewrite the metric of the 
spanned by r and the as 

dr"^ + r'^dnl = dp^ + p^d^l + dy^ + y^dcj)^ , (2.7) 

where = cos^ 9 and = sin^ 6 and hence r'^ = + y"^. The D7-branc is extended 
along p and the S^. We have used polar coordinates y and <p for the {xs,xg) plane. On 
the D7-brane worldvolume y and appear as scalar fields, which generically depend on 
all the D7-brane coordinates (^^, . . . ,^^). The mode of y that is independent of the 
coordinates is dual to Om- To describe a constant real mass M, preserving space and time 
translation symmetries, we need a solution of the form y{p). The solution to the equations 
of motion is then simply a constant, y(p) = L, which represents flavor fields with a mass 
M = L/(27ra') [26,29]. 

In the trivial case, y{p) = 0, the D7-brane clearly preserves the SO{4) x ^7(1)^ isometry 
of and (j), which is dual to the symmetry preserved by a massless hypermultiplet. In 
the nontrivial case, y{p) = L, the D7-brane preserves only the iS'0(4) isometry, dual to the 
symmetry preserved by a massive hypermultiplet. 

The trivial solution describes a D7-brane that fills AdS^ x S^. Massless flavors in the 
probe limit will indeed preserve conformal invariance, hence the appearance of the AdS^ 
factor. Recalling that the AdS5 radial coordinate is dual to the CFT energy scale, we 
can interpret the non-trivial solution y{p) = L as the renormalization-group (RG) flow 
associated with massive flavors. Here the induced D7-brane metric is AdS^ x S*^ near the 
boundary p — >■ oo, but as the D7-brane extends into AdS^, the shrinks and eventually 
collapses to zero size (smoothly) at p = 0. Recalling that the AdS5 radial coordinate is not 
p but r = \J fP" + y{pY, we see that at p = the D7-brane has only reached r = L. Prom 
the perspective of an observer in AdS^, the D7-brane simply ends there. In the RG flow, 
the flavor flelds are absent at energy scales below their mass gap. 

The spectrum of fluctuations about the non-trivial solution y{p) = L was computed 
in ref. [29]. In that case the fluctuations are normal modes, i.e. standing waves trapped 
between the AdS^ boundary and the endpoint of the D7-brane. These normal modes are 
dual to poles in the two-point functions of the dual operators, i.e. mesons, which in the 
large- A'^^: and probe limits are exactly stable. These mesons do not arise due to chiral sym- 
metry breaking or conflnement, rather they arc simply deeply bound states, with masses 
on the order of M/\/A [29]. We will not review the meson spectrum in detail here, but 
we will mention that the mesons form representations of the 50(4) symmetry, and in fact 
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because the mesons form short multiplets of the J\f = 2 superconformal algebra, their di- 
mensions are completely determined by their R-symmetry charges [29] . The dual statement 
is that fluctuations of D7-brane fields can be decomposed into spherical harmonics of the 
S^, which cannot mix since they form different representations of SO{A), and higher 
angular momentum implies larger (Kaluza-Klein) mass, which means the dual operator has 
a larger dimension. 

2.1 Finite Density States 

The conserved current associated with baryon number U{1)b, J", is dual to the compo- 
nents of the worldvolume gauge field in the Minkowski directions that are independent of 
the coordinates. Ay [27]. In particular, if we want to study states with a finite chem- 
ical potential ^ and charge density (J*), preserving space and time translations, we must 
introduce a worldvolume At{p). We can also motivate the introduction of At{p) with the 
following physical argument. A single charged excitation is represented in the bulk by a 
single string that ends on the D7-brane. To obtain a charge density, we must introduce 
a density of strings, which will source the worldvolume gauge field, producing an electric 
field in the radial direction, Fpt{p). If we choose Ap = Q gauge, then to obtain Fpt{p) we 
must introduce At{p). 

Notice that a number on the order of N^. strings is strong enough to pull the D7-branes 
all the way to the "bottom" of AdS^^ meaning r = 0, hence any solution representing a 
finite baryon density of order^ Nc will not end, i.e. the will not collapse at finite r [27]. 
In such cases the fluctuations of worldvolume fields will not be normal modes but quasi- 
normal modes (QNM's), that is, the eigenfrequencies will become complex, since some part 
of the standing wave can now "leak" across the (Poincare) horizon of AdS^. That conforms 
to our field theory expectation that any finite density will give mesons a finite width. 

Luckily, two exact solutions describing finite M and p are known [28]. The first is 
simply At{p) = p with constant y{p) = L. These constant fields are a solution of the 
equation of motion because the field strength Fpt = 0, so the gauge field equation is 
trivially satisfied and the scalar equation is unchanged from the At{p) = case. These 
solutions describe zero density states, however, and so are only relevant when studying the 
grand canonical ensemble. The second solution can be written in terms of an incomplete 
Beta function, 




where Af = Tmlir"^ = , and the parameters c and d are fixed to be 

c = -fJ\fi2TTaf {p^ - M^) M, (2.9) 

d = 7AA(27ra')4 {p'^ - M^) p, 

^Following ref. [28], we employ a normalization in which a single charge/string has baryon number one, 
so that a baryon actually has baryon number N^. 
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with 7 = [^B -j)] ^ « 0.363. These are related to the expectation values of the dual 
operators via (Om) = {2iTa')c and (J*) = d. We have also introduced the parameter 

c(2W)_(Oa/)_M 



Notice that y{p = 0) = 0, indicating that the S"^ collapse occurs only when r = \/p^\-y{p)^ 
is zero, that is, at the bottom of AdS^. Notice also that although M = implies c = and 
hence y{p) = 0, the gauge field At{p) = \y{p) remains nontrivial. 

For equilibrium states, the AdS/CFT dictionary equates the on-shell bulk action with 
the field theory grand canonical partition function. By comparing the free energies associ- 
ated with the two solutions above, we discover a phase transition [28]. When p < M, the 
thermodynamically preferred solution is the trivial one, describing states with {Om) = 
and zero density, (J*) = 0. When p > M, the preferred solution is the non-trivial one, 



describing states with the nonzero {Om) and (J*) given by the c and d in eq. (2^). The 
critical exponents associated with the transition take mean-field values. Technically, this 
zero-temperature transition is a quantum phase transition, albeit a rather boring one from 
nothing (zero density) to something (finite density). Notice that the transition is only 
apparent in the grand canonical ensemble. In the canonical ensemble we have access only 



to finite-density states, by definition. Notice also that the parameter e defined in eq. ( 2.1C| ) 



is bounded: < e < 1. The limit e = corresponds to zero mass, or finite mass and 
infinite chemical potential. The limit e = 1 corresponds to M = ^u, where the transition to 
zero-density states occurs. 

Upon introducing a finite temperature, dual holographically to studying the same 
system but in an AdS-Schwarzschild black hole spacetime, straightforward calculations 
reveal two unusual features at temperatures low compared to the density. First, the heat 
capacity scales with temperature as cy oc T^, in stark contrast to a Fermi liquid, where 
cy oc T, or a Bose liquid, where cy oc in three spatial dimensions [30]. Second, the 
zero-temperature limit of the entropy density s is nonzero: 

limr^os = ^\/A(J*) [30]. As 

argued in ref. [32], these are essentially the heat capacity and entropy associated with 
a single charge/string, times the number of charges /strings, (J*). The low-temperature, 
finite-density states appear to be locally thermodynamically stable [52]. 

To our knowledge, the only analyses of the fluctuation spectrum about the zero- 
temperature, finite-density solution appear in refs. [30,31,53],^ where a pole was discovered 
in the two-point function of the U{1)b charge density (and current) two-point functions, 
which for massless flavors M = has a dispersion relation^ 

1 

uo(k) = ±—k + i— + O (k^) , (2.11) 
V3 6/i 



^With finite baryon density and finite temperature, various parts of tlie fluctuation spectrum have been 
computed in refs. [54-62]. 

*Our frequency and momenta, ui and k, have signs opposite to those of ref. [30]. In our conventions, the 
zero sound appears in the upper-half of the complex uj plane, hence the imaginary part has a sign opposite 
to that of ref. [30]. 
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which is formally identical to that of Landau's "zero sound" pole in the density two-point 
function in a Fermi liquid, hence the mode was dubbed zero sound by analogy. Notice in 
particular that the speed is identical to the speed of sound in a finite-temperature conformal 
plasma in (3+1) dimensions. 

Our objective is to analyze fluctuations of all bosonic worldvolume fields, including 
finite-momentum fluctuations, in order to detect perturbative instabilities. Before pro- 
ceeding to the details of the calculation, let us review the circumstantial evidence in favor 
of instabilities. 



2.2 Reasons to Expect Instabilities 

We have a number of reasons to expect instabilities, based on intuition from similar theories 
at weak coupling and also from similar holographic systems. 

The first reason is simply the presence of a nonzero entropy density at zero temperature, 
which renders the system vulnerable to instabilities. 

The second reason is the presence of charged fermions in our field theory. Our intuition 
here comes from QCD with an asymptotically large baryon chemical potential. In such a 
regime we expect Fermi surfaces of weakly-interacting quarks (technically Landau quasi- 
particles). When Nc = 3, standard arguments from the BCS theory of superconductivity 
suggest that the quarks form Cooper pairs, producing a color superconductor.^ In the 
large-A'^c limit, however, things may change. The quarks still interact weakly, but the 
diagram describing the relevant attractive interaction between them, providing the "glue" 
for the Cooper pairs, is non-planar, and the superconducting gap goes as e"^'^/^, so both 
are suppressed in the 't Hooft limit [36,37]. One proposal for the actual ground state [36] 
is for quarks to pair with holes in the Fermi surface, moving in parallel with momenta 
on the order of fi, producing a gauge-invariant condensate, the usual chiral condensate, 
but spatially modulated with a wave vector of magnitude 2fi. In our theory, the M = 2 
hypermultiplet includes fermions charged under U{1)b which may be susceptible to such an 
instability. Whether our finite-density states include a Fermi surface is unclear, however: 
the system has a zero sound mode, like a Fermi liquid, but a heat capacity unlike a Fermi 
or Bose liquid. 

Inhomogenous ground states have also been found in several models of holographic 
matter, including probe brane systems [38-43]. In these cases, the crucial ingredient was a 
Chern-Simons term in a (4+l)-dimensional bulk spacetime [38-41] or an axion-like coupling 
in a (3+l)-dimensional bulk spacetime [42,43]. Our model includes Chern-Simons-like 
terms via the D7-brane WZ terms, providing yet another suggestion that our system may 
have a finite-momentum instability. We should mention, however, that in our analysis 
below the WZ terms only enter the equations of motion for the D7-brane gauge field 
components in the directions, not in the AdS^ directions (see eq. ( 3.6fD ). 

Another major reason to expect instabilities is the presence of charged scalars in our 
field theory, for which a chemical potential generically triggers Bose-Einstein condensation. 
The canonical example of Bose-Einstein condensation is a massive complex scalar, which 



^For a review of the general arguments, see ref. [63] and references therein. 
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necessarily has a global U{1) symmetry, with quartic self-interactions. A chemical potential 
for the U{1) effectively acts as a negative mass-squared. A chemical potential larger than 
the mass triggers a second-order phase transition: the potential develops a new minimum, 
the scalar acquires an expectation value, and the U{1) symmetry is broken. 

In our field theory, the hypermultiplet contains two complex scalars, the squarks, which 
indeed have quartic self- interactions, as well as Yukawa and other interactions [49]. The two 
squarks q and are in the Nc of SU{Nc), form a doublet of SU{2)ji, carry identical U{1)b 
charges, and are neutral under SU{2)l x If the squarks acquire an expectation 

value, then we expect SU{Nc) x SU{2)ji x U{1)b to be broken to some subgroup. 

What would squark condensation look like holographically? Consider the intersecting 
branes in table (|2.2| ). Upon taking the usual Maldacena limits, we obtain AdS^ x S^, 
where we may imagine that the D3-branes are "hidden" behind the AdS^s Poincare horizon. 
We then introduce Nf D7-branes with worldvolume U{1) electric flux. The D3-branes 
feel two forces in the radial direction, one towards the horizon, due to gravity (AdS is a 
potential well), and one towards the boundary, due to the electric flux. If the electric flux 
is sufficiently large, then the D7-branes may "suck" D3-branes out of the stack of N^.- The 
authors of ref. [64] called this "Fermi seasickness." 

In our case, the D3-branes may then dissolve into the D7-branes, becoming instantons. 
The instanton moduli, namely the position, orientation in U{Nf), and size, are isomorphic 
to the scalar expectation values, including those of the squarks and the adjoint scalars of 
= 4 SYM [65-68]. Specifically, the size modulus is isomorphic to the squark expectation 
value. Normally, for a single supersymmetric D7-brane, the U{1) worldvolume theory 
admits only point-like instantons, that is, instantons with zero size modulus, dual to squarks 
with zero expectation value. In the presence of nonzero background electric fiux, however, 
the instanton could possibly acquire a finite size.^*^ The instanton number counts the 
number of D3-branes that have dissolved, and hence indicates breaking of SU{Nc)- a 
single-instanton solution implies the breaking SU{Nc) — )• SU{Nc — 1) x U{1). We thus 
expect a squark condensation instability to appear in the bulk as the D7-brane sucking 
D3-branes out of the stack, which then dissolve into the D7-brane, becoming finite-size 
instantons. Crucially, however, the rate at which the D7-branes suck D3-branes out of the 
stack goes to zero when Nc ^ oo and Nj <^ Nc [64], that is, the process takes forever. We 
will thus not be able to detect this instability, since we work in the probe limit. 

What symmetry-breaking instabilities can we detect, in principle, from our holographic 
calculation? We focus only on fiuctuations of D7-brane worldvolume fields, dual to mesonic 
operators. These are invariant under SU{Nc) and U{1)b, so we cannot detect breaking 
of these symmetries. The various mesons of our theory form various representations of 
the remaining global SU{2)l x SU{2)b symmetry, and the U{1)r symmetry present when 
M = 0, so we may detect instabilities toward states in which these symmetries are broken 

^"The possibility of finite-size instantons when Nf = 1 should generically occur whenever the two stacks 
of branes feel a force towards one another. For example, a nonzero Kalb-Ramond B-field, which introduces 
non-commutativity on the worldvolume of the D7-branes, will also attract the D3-branes to the D7-branes, 
and allows for finite-size instantons even when Nf = 1. In the holographic context, the B-field is dual to a 
Fayet-Iliopoulos term in the field theory [69]. We thank D. Tong for discussions on this point. 
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to some subgroups. We may also detect instabilities towards states breaking translational 
or rotational symmetries, for example if we see an instability towards condensation of vector 
mesons, or a finite- momentum instability. 

In AdS/CFT we only have access to gauge-invariant observables. If we do see an 
instability, then most likely we will not be able to identify the microscopic mechanism. 
For example, in principle we could detect an instability towards a state in which {Om) 
has a plane- wave form with wave vector of magnitude 2|U, just like the chiral density wave 
in large-A'^c QCD. We would probably not be able to determine whether the microscopic 
origin of the instability in the field theory was the same as in large-A'c QCD, that is, 
whether quarks were pairing with holes in a Fermi surface. AdS/CFT provides access only 
to "macroscopic" data, namely correlation functions of gauge-invariant operators. 

Nevertheless, a microscopic interpretation of the solution in eq. ( p.^ ) was proposed in 
refs. [28,48]. The suggestion is that the solution in eq. (|2.8| ) in fact represents squarks 
that have already condensed. The heuristic field theory argument goes as follows. We 
begin with nonzero M but /i = 0, where the density (J*) = and also {Om) = 0. We 
then increase ^ until /i > M, at which point the scalars condense, producing a nonzero 
density (J*) and nonzero {Om)- Indeed, the form of {Om) in eq. ( |2.9| ) is identical to that 
of a charged scalar with quartic coupling. Moreover, we may anticipate that the system is 
perturbatively stable because, in some sense, the instability already took place: the squarks 
Bose-Einstein condensed. 

A number of questions arise about such an interpretation, however. As mentioned 
above, if we trust our intuition from the canonical example of Bose-Einstein condensation, 
then if the squarks condense we expect SU{Nc) x SU{2)fi x U{1)b to be broken to some 
subgroup. The solution in eq. ( p.^ ) describes a state in which the only operators that 
acquire expectation values, J* and Om, are neutral under these symmetries. The dispersion 
relation in eq. ( |2.11| ) looks suspiciously like that of a Goldstone boson, associated with 
U{1)b breaking for example, but as argued in ref. [21], such a Goldstone boson would have 
a damping coefficient suppressed by a factor of Nc relative to that of eq. ( p.ll| ). In short, 
the solution in eq. ( |2.8| ) preserves all global symmetries. The solution of eq. ( p.8| ) also 
has instanton number zero, that is, the D7-brane has zero D3-brane charge, and hence no 
instanton size modulus. That suggests the squarks have zero expectation value and that 
SU{Nc) is not broken. How can such a solution represent squark condensation? 

Of course, our intuition may fail, due to an essential difference between our theory and 
the canonical example of Bose-Einstein condensation: the squarks are not gauge-invariant, 
so if they develop an expectation value, identifying the broken symmetries requires some 
caution. The SU{Nc) gauge symmetry of our theory is not a physical symmetry at all, but 
merely a redundancy of our description, so no physical, gauge-invariant order parameter 
exists that can detect its breaking. In a weakly-coupled theory, such as the electroweak 
theory, a scalar expectation value merely serves as a "convenient fiction" to identify the 
breaking (to quote ref. [63]). A physical statement of spontaneous breaking of gauge 
invariance is, for example, the presence in the spectrum of a longitudinal mode for a vector 
boson. To detect breaking of a global symmetry, we must find a gauge-invariant charged 
operator and demonstrate that its expectation value is nonzero. For example, we can 
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use the scalar expectation value itself to construct a gauge-invariant charged operator, by 
contracting its gauge indices with copies of itself [63] . 

We may thus be able to argue that, in fact, the squarks can condense without breaking 
any symmetries. Consider for example the U{1)b symmetry. In our theory the only gauge- 
invariant operator charged under U{1)b is the baryon operator itself,^^ which due to anti- 
symmetry in color indices is actually constructed from Nc quarks, hence we cannot build 
a baryon from a squark expectation value. We thus learn a very important lesson: in 
our theory, a nonzero squark expectation value does not immediately imply that U{1)b is 
broken, a fact that conflicts with our intuition from the canonical example. To be clear, 
our point is not that U{1)b can never break, but rather that a nonzero squark expectation 
value does not provide sufficient information to demonstrate U{1)b breaking. 

Even if we could somehow argue that the squarks condense without breaking any 
symmetries, we should ask how a system of condensed scalars can have a zero sound 
excitation and a heat capacity cy cx T^. Both of these features seem strange, if we trust 
our intuition for weakly-coupled scalars. Of course, again our intuition may fail: the zero 
sound and unusual heat capacity may be strong-coupling effects. 

Determining what microscopic state the solution in eq. (|2.8| ) actually represents is 
beyond the scope of this paper, but we can further characterize that state by asking about 
the spectrum of fluctuations. We now proceed to the calculation of that spectrum. 

3. Fluctuations of D7-brane Fields 

We want to calculate the spectrum of fluctuations of bosonic D7-brane worldvolume fields. 
To do so we must determine the equations of motion for the fluctuations about the back- 
ground solution described in section Our analysis will closely follow that of ref. [29], 
where the spectrum of fluctuations about the constant solution y{p) = L with At{p) = 
was computed, corresponding to the meson spectrum for finite-mass flavors in states with 
zero density. 

We will employ a very general ansatz for the fluctuations, allowing them to depend 
on the Minkowski directions as well as the radial direction and the angles on the S^. 
Specifically, we introduce a fluctuation of the embedding, y{p) — >■ y{p) + (27ra')x, of the 
other worldvolume scalar, (p = {2-Ka')ip, and of the worldvolume gauge field. A), = d^At + 
5A}y. The corresponding fluctuations of the induced metric and the U{1) gauge field, up to 

Actually, the baryon is the only IocmI operator charged under U{1)b- Non-local operators charged under 
U{1)b should also exist. We thank A. Cherman for discussions on this point. 

^^If Nf > Nc, we can construct a baryon using the squarks by anti-symmetrizing in both color and flavor 
indices. We have Nf = 1 <^ Nc, however. 

^^A calculation of the back-reaction of D7-branes, with (7(1) electric flux, on the flelds of supergravity 
appears in ref. [70]. Although apparently the ansatz used in that calculation was incomplete, as shown in 
ref. [71], the results suggested that the electric flux produces a potential for the size modulus of D7-brane 
instantons, in keeping with the fleld theory expectation that the U{1)b density should produce a nonzero 
squark expectation value. The authors of ref. [70] argued, however, that U{1)b is not broken when the 
squarks acquire an expectation value, which is consistent with our arguments. 
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quadratic order, take the following form: 



ab 



s 



G 



fii/ 



6Fab = da5Ab - db5Aa. (3.1) 
We will only need the the pullback of the four- form to first order in the fluctuations, 

5P[C^^'^] = cos^ e{27ra') (^dp (^^^ dp A dfls + ^dt^dt A dfl^ + ^V^dx A d^s^ . (3.2) 



Inserting eqs. (^) and (3^) into the D7-brane action eq. ( |2.6| ), we find that the 
terms linear in fluctuations vanish because the background fields are on-shell, as expected. 
To write the terms quadratic in fluctuations, S'j^j^ \ let us first define, in terms of the 
background solution. 



/_ ^2 



Vab 



gab + (27ra')-Fa 



e 



al3 



\ 



Tgij/ 



(3.3) 



where a, (3,^ denote Minkowski directions, i,j,k denote 5^ directions, with gij the metric 
of the 5'^, and we have ordered the coordinates as {t,p,a,i). The inverse of rjab will be 
denoted r]"'', for which we will define symmetric and anti-symmetric parts as 



(3.4) 



Notice that the only non- vanishing element of rj'^ is ry^* = — The terms in the action 
quadratic in fluctuations are then 
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^ .ba dc ^ da „bc 



T] I —rj rj 



-r/ ?? ) X 



^{daxS'b + dbx6P) + da5Ab - dbSA^ 



+r\dp6A, - didAp)djdAke'^'' 



(27ra') 



MdcXS'a + ddxS'c) + dcSAd - d^SA, 



cos^ e^da^idpdA^ - d^6Afs)e''^"', (3.5) 



where we have dropped overall constant factors, including for example a factor of {2iTa')'^T£)j, 
since these will not affect the equations of motion. Varying eq. ( |3.5| ), and then making the 
gauge choice 5Ap = 0, we find the equations of motion for 99, X) ^t, Aa, and Ai, 
respectively: 



1 



daiV^Vs^—^db^) = 0, 

pZ _|_ yZ 



(3.6a) 
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[p^ + y 



2\2 



p2 _|_ y2 



_|_ yZ 

-^dpdtx + da{^r]PpT]fdpbA^) = 0, 



p2 _|_ y2 



dhX) 
(3.6b) 

(3.6c) 



_|_ yA 



-da{ 



^rf'rifdtSA,) 



0, 



pA _|_ yi 



-7]7]''r^fdb6At) 



(3.6d) 



0(3.6e) 



y 



P^ + y 



:dtdjx - da{^/^Vsr]''dj6Ah) + da{V^v''VsdbSAj 



-4ip^+y^){p + yyy^^dj6Ak = 0. 



(3.6f) 



With our choice of gauge, the equation of motion for 5Ap, eq. ( p.6c| ), becomes a constraint 
on the remaining fields. Notice that the final term in eq. ( |3.5| ) does not contribute to the 
equations of motion, due to the anti-symmetry of e°^^'^. In other words, the WZ terms 
only affect the gauge field fiuctuations Ai. That is in contrast to ref. [38], for example, 
where the Chern-Simons terms entered the AdS^ gauge field equations of motion. Notice 
also that if we set Fab = in our r/^fe and also set y'{p) = 0, then our equations reduce to 
those of ref. [29], up to a difference in gauge: we use 5Ap = gauge, while in ref. [29] the 
gauge choice —dtdAt + daSA^ = was used. 

Given that the background solution preserves translational and rotational invariance 
in the Minkowski directions, we will assume without loss of generality that the fiuctuations 
have a plane- wave form with momentum only in the direction^^, e*'^*-«'=^3_ jj^ that case, 
the transverse gauge field fiuctuations 5Ai and 6A2 decouple from the other fluctuations, 
and in fact obey exactly the same equation as if, eq. (3.6a). The remaining fiuctuations 
remained coupled to one another. 

Following ref. [29], we will decompose the fiuctuations into spherical harmonics and 



look for ansatze that decouple fields from one another. Let and 3^, 
vector spherical harmonics on S^, which obey the following relations: 



m,± 



m,± 



-m(m + 2)y™, 
-(m+l)2y-'±. 



be scalar and 



(3.7) 



m.it 



0, 



^^Notice our choice of signs: we use e'"' ''^^^ rather than e »"'+*'=^3^ With our conventions, stable modes 
wiU appear in the upper-half of the complex cj plane. 
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where V*, i?*^ = 25*^, and g are the Laplacian, Ricci tensor, and determinant of the 
metric. With respect to the SO{A) = SU{2)ji x SU{2)l isometry of the 5^, the scalar 
harmonics 3^™" transform in the (^j representation, with m > 0. We can construct a 
vector harmonic from the scalar harmonic by taking VjJ^™'. The vector harmonics y^' 
transform in the (^^^^, ^^^^) representation, with m > 1. We can decompose the scalars 
and X) and the gauge field components A^, into S'^ scalar spherical harmonics, and the 
gauge field components Ai into vector spherical harmonics. 

The ansatze that ultimately decouple the fluctuations are 



k 



^iuit—ikx^ -ym 

X = 6At = 5A^ = 0, 6Ai = $±(p)e^'^^-*^^3y 
X = 5At = 0, M3 = $(/5)e^'^*-^^'^33^"', 5A 



^iuit—ikx-j -V im,it 



(3.8a) 
(3.8b) 

^(^py^t-ikx^y^ym^ (3.8c) 



where rj^^ is defined by rj^^ = rj^^g^^ = ^g '^K The ansatz in eq. ( 3.8a| ) involves only a 
single function <&(p). The ansatz in eq. ( 3.8b| ) involves only fluctuations with purely vector 
harmonics on the , which decouple from fluctuations involving scalar harmonics since 
these are in different representations of 50(4). In contrast, the ansatz in eq. ( ^.8c| ) involves 
a coupling between 5^3, a scalar harmonic on the , and the vector harmonics on the 
built from VjJ^"^. The ansatze for 5Ai, and 5A2 are the same in form to that of x iii 
eq. ( p.8a ) , that is, ip = $(p)e*'^*~*'^^f3^™', and similarly for 5Ai and 5A2. 

Inserting these ansatze into the equations of motion, eqs. ( 3.6a| ) to ( |3.6^ ), and using 



4(> 



— ab 

-ms 2 



1 



;dbx) = 0, 



(3.9a) 



w'=S*'5,Mz) - daiV^v''VsdcSA,) + 4(p2 + y^)(p + yyy^''dj5Ak = 0, (3.9b) 
^dbSAs) = 0. (3.9c) 



Eqs. (|3^ ), (l3!9b|) , and (|3^) are the non-trivial equations resulting from the ansatze in 
eqs. (|3.8a| ), (|3.8bD , and ( ^.8c ), respectively. Eqs. ( 3.9^ ) and ( 3.9c ) are identical in form, 
and indeed are identical in form to the equations of motion of ip, 5Ai, and 6A2. We have 
thus found that both of the scalar fluctuations and the fluctuations of the gauge field in 
Minkowski directions all obey the same equation of motion. Notice that these fluctuations 

built from y"" and/or V^y". The fluc- 



i.e. 



are built from scalar harmonics on the S 
tuations built from genuine vector spherical harmonics 3^™'=*= obey a different equation of 
motion, eq. ( |3.9b ), which will actually produce two equations, for ^^{p) and <I>~(p). 

We thus have three "master equations" for the fluctuations. Let ^{p) collectively 
denote the /?-dependent part of the two scalar and the Minkowski gauge field fluctuations, 
and ^^{p) denote the p-dependent part of the vector harmonic fluctuations. Using the 
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relations for the spherical harmonics, we can write the master equations as 

dpiV^'qPPrf'-'dp^) - ^ (tj2^"r/^^ + ^^(r,^^)^ + m{m + 2)7/^=^7?^^) $ = 0, (3.10a) 

±A{p' + y^){p + yy')(m + l)^) $± = 0, (3.10b) 



where the first equation comes from eq. (3.9a) and the second comes from eq. ( 3.9b| ). 

Recalling that the AdS^ radial coordinate corresponds to the field theory energy scale, 
we may interpret the background solutions y{p) and At{p) as an RG flow. The fact that 
some of our fluctuations couple corresponds to the fact that the dual operators mix under 



that RG flow. For example, in eq. ( 3.8a ) we see that Xi and 5A^ couple to each other, 
but decouple from all other fluctuations, hence in the field theory the dual operators Om-, 
J*, and J^^ mix with each other, but do not mix with any other meson operators. 

To show that we have included all fluctuations, let us compare to the zero-density 
accounting of fluctuations of ref. [29]. At zero density, both of the scalars decouple from 



all other fluctuations. As is clear from eq. (3.6a), our (p decouples. For x things are more 
subtle. Suppose we set the mass M = 0. In that case e = 0, and from eq. ( p. 8a ) we see 
that indeed x decouples from all other fluctuations. The reason is that when M = the 
chiral f7(l)_R symmetry is restored, and (p and x are the only worldvolume fields charged 
under that symmetry [29]. In other words, they are in a different representation of f7(l)_R 
from all the other fluctuations, and hence must decouple. When M is nonzero, the U{l)ii 
is broken. Nothing then forbids x from coupling to other fluctuations, so we must allow 
for fluctuations of the form in eq. (|3.8aD, although in the end x does decouple, as we see 



in eq. ( 3.9a ). The upshot is that eq. ( 3.8a ) is the generalization to finite density of the x 
fluctuation in ref. [29]. 

In ref. [29], three types of gauge field fluctuations were identified, which were called 
type I, II, and III. Type I were identical to our eq. ( |3.8b ), and our eq. ( |3.8c ) is essentially the 
generalization of the type III fluctuations to flnite density and to our choice of gauge. The 
type II fluctuations involved scalar harmonics of the Minkowski components of the gauge 
fleld, which decoupled from all other fluctuations. In our case, where the charge density 
breaks Lorentz symmetry to spatial rotations and translations, the Minkowski components 
of the gauge field split: 5Ai and 6A2 still decouple, but as we discussed above, when M is 
nonzero 5 At and SA-^ mix with x- 

Notice that m = is a special case. When m = 0, the ansatz in eq. ( ^ ) is ill-defined. 
That is not unexpected, since the type III fiuctuations of ref. [29] were also singular when 
m = 0. For the ansatz in eq. ( |3.8a| ), eqs. ( |3.6b ) and ( p. 611) are constraint equations, 



coming from setting 5Ap = and 5Ai = 0, respectively. When m = 0, eq. ( p.6fD is trivially 
satisfied, meaning one constraint is eliminated. We then expect to find an additional degree 
of freedom, which in our case is the fiuctuation giving rise to the zero sound mode, as we 
review in the appendix. 

We will make two final simplifications, to make an analysis of the singular points of 
the equations of motion easier, and also to make our numerical analysis easier. First, we 
will define dimensionless variables. Let N = Af/c, which has length dimension three (recall 
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eq. (p.9D ), which we use to define the dimensionless variables /?, y, u, and k: 

ff = P'N' -iy\ y'^ y'N' -iy\ (3.11) 
-^--'(^v(^-i))''', Ai^ = (ivV - ''''' 

In terms of these new variables, we have 



1/1 ,Y^^^ „ f 1 lA dy 



6Ve2 y Vp' + 1 6 3y dp ^(^6 + 1) (^ - 1) 

The three master equations now take the form 

^" +Pi{pW + qi{p,oj,k,m)^ = 0, (3.13a) 
+P2{p)^^' + qt{p,0,k,m)^^ = 0, (3.13b) 



where 



Flip) = (3.14a) 
1 + p^ 

. r ^ -m(m + 2)/>4 ^^(1+^') 

'^^^^''^'^'"'^ = i + p-6 + (p2 + y-2)2 - JYTWWTW' ^ ^ 

p^(2 + 5p6) + (-2 + f)f + 4p(l + p6)yy^ 
P2 p) = -n ^ -&\(-2 I -2^ ' 3.14c 



(p2+y2)2 (l + p6)(^2+y2)2- 

The master equations, eqs. (|3.13a|) and (l3.13bD , and their solutions, will be the focus of 
the remainder of this paper. 

Finally, let us demonstrate how the fluctuations "see" a near-horizon AdS2 region 
[21,33-35]. In our case "near-horizon" means p — )• 0. In that limit we have 

m = ^^+0{f). (3.15) 



Consider the equation for $(p), eq. ( p. 13a ). From eq. ( p. 14b ) we see that as p — t- 0, the 
dominant term in the coefficient qi{p,Ljj,k,m) is the one proportional to u)^, which will 
scale at leading order as 1/p^, while the term proportional to m{m + 2) scales as p'^ and 
the term proportional to k'^ scales as p^. Here we already see a hint of AdS2- only the 
radial and time directions will appear in AdS2, so we should see a suppression of terms 
involving spatial momenta. Here we actually see a suppression of terms involving the linear 
momentum in Minkowski directions as well as the angular momentum in the 5^ directions. 
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Prom eq. ( p.l4a| ) we can see that when p — t- the coefficient pi{p) — )■ 3p^, which we will 
treat as sub-leading, i.e. we will drop the term from eq. ( |3.13a| ) in the near-horizon 
limit. The equation for $(/o) in the /) — t- limit thus becomes, at leading order, 

$" + ^^(l-e2)$ = o. (3.16) 

If we extract a factor of p, that is, if we define <!>(/)) = p^{p), then eq. (|3.16| ) becomes 



-,2 

-2\ 



2 id 

|," + _|,' + _(l_e2^$ = 0, (3.17) 

which is precisely the equation of motion for a massless scalar in AdS2, with metric ^ — 
p^dt^ , and with frequency a)\/l — e^- 

The p — )• limit for the 'I*^(p) equation, eq. ( p.l3b| ), is even simpler. From eq. ( 3.14d| ) 



we see that as p — t- the leading term in q2{p,u}, k,m) is identical to the leading term in 
qi{p,Lo,k,m), namely tD^(l — e^)/p^. From eq. (|3.14c ) we find that p2(p) has an especially 



simple behavior as p — )• 0, namely P2{p) — ^ 2/p. The equation for the *I'^(p) in the p — )• 
limit thus becomes, to leading order, 

p p4 

which again is precisely the equation of motion for a massless scalar in AdS2 with frequency 
a)\/l — Here we did not even need to rescale $^(p) by a factor of p, in contrast to $(p). 

We emphasize that the appearance of a near-horizon AdS2 region is due to the presence 
of a nontrivial At{p), or in field theory language a finite density. At zero density, Lorentz 
invariance demands that the equations of motion depend only on the Lorentz-invariant 
quantity uj'^ — k'^, hence dependence on k could not be suppressed relative to w in the 
near-horizon limit. 

Given that eqs. ( p. 17 ) and ( p.l8| ) describe a massless scalar field in AdS2, for all 



values of k, e, and m, we will not see an instability involving a violation of the AdS2 BF 
bound, that is, an instability arising from the mass-squared of a field in AdS2 becoming 
too negative. Of course, many other kinds of instabilities are possible. In fact, for our 
system two deformations are known to produce instabilities without any violation of the 
AdS2 BF bound. The first is a sufficiently large U{1)b magnetic field, which for massless 
flavors causes spontaneous breaking of the chiral U{1)r symmetry [48,72,73]. The second 
is an isospin chemical potential, meaning we introduce a second D7-brane coincident with 
the first, so that the field theory flavor symmetry is enhanced to U{1)b x SU{2), and we 
then introduce a chemical potential associated with the Cartan generator of SU{2). That 
causes vector meson condensation, leading to a p-wave superfluid state [74-77]. In each 
case a near-horizon AdS2 is present, but no violation of the AdS2 BF bound occurs. 

4. The Spectrum of Quasi-Normal Modes 

Every fluctuation of a D7-brane worldvolume field is dual to a meson operator. The 
AdS/CFT dictionary equates the on-shell bulk action with the field theory generating 
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functional. The solutions of the bulk linearized equations of motion, eqs. ( |3.13a ) and 
(3.13b) thus determine the dual two-point functions. In Lorentzian signature we have a 
choice of two-point functions: retarded, advanced, etc. In the bulk these are determined 
by the boundary condition we impose in the near-horizon region. We want the retarded 
two-point functions, so we will impose in-going boundary conditions in the near- horizon 
region, i.e. the solutions will look like waves traveling towards the bottom of AdS^ [78]. 

We will not compute the full two-point functions. To study stability, simply computing 
the eigenfrequencies of the linear operators in eqs. (|3T3a| ) and ( |3T3b| ), that is, the QNM's 
of the system, will be sufficient. These QNM's are dual to poles in the retarded Green's 
functions. An instability would appear as a pole in the lower half of the complex frequency 
plane^^, since such a mode would have uncontrolled exponential growth. 

The main technical difficulty in computing the QNM's is the presence of irregular 
singular points in eqs. ( 3.13a ) and ( 3.13b| ) . To deal with these, we will employ a combination 
of analytical and numerical techniques. In section 4.1 we consider the low-frequency limit, 
where we can match asymptotic expansions to determine analytically the low-frequency 
form of the two-point function, following ref. [16]. In section 4.2 we discuss the singular 
points of eqs. ( |3.13a| ) and (3.13b), review Leaver's method for dealing with these, explain 
why Leaver's method is difficult to use in our case, and then present our own numerical 



method, which we use to compute the QNM spectra presented in section 4.3 



4.1 Lovif- Frequency Expansion 

We expect an instability to appear as a QNM crossing into the lower-half of the complex 
frequency plane at small frequency. Luckily, we can obtain the exact form of the retarded 
Green's functions in the low- frequency regime, following ref. [16]. We divide the radial 
direction into two regions, the "inner" and "outer" regions, defined by p ^ 1 and 0,/p <^ 1, 
respectively. These regions overlap at small frequency Q <^ 1. The inner region is simply 
the near-horizon AdS2 region. 

The presence of the near-horizon AdS2 suggests that in the field theory in these finite- 
density states some (0 + l)-dimensional CFT emerges at low energies. We may invoke 
an AdS/CFT correspondence for the AdS2 region itself: every field in AdS2 will be dual 
to some operator in the dual (0 + l)-dimensional CFT, living at the boundary of AdS2, 
where we match to the outer region. For example, if a bulk field in the AdS2 region obeys 
the equation of a massless scalar, then that field is dual to some marginal (dimension 
one) scalar operator in the (0 -|- l)-dimensional CFT. We expect the Green's functions of 
the (0 + l)-dimensional CFT to determine the low-frequency behavior of the full Green's 
functions [16], which indeed occurs, as we will show. 

In the last section we saw that all bosonic fiuctuations of D7-brane fields ultimately 
obey the same equation of motion in the AdS2 region, that for a massless scalar, regardless 
of their spin (scalar or vector), their linear momentum in Minkowski directions, or their 
angular momentum in directions. We thus conclude that in the field theory, every 

Recall again our sign convention, such that stable modes appear in the upper-half of the complex 
frequency plane. 
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bosonic meson operator is dual at low energies to some marginal scalar operator in the 
(0 + l)-dimensional CFT, regardless of its spin, momentum, or its SO{4) charges. 

We begin with $(/o). The equation of motion in the inner region is eq. ( ^.16 ). Extract- 



ing a factor of p via ^{p) = p^ip), we obtain eq. ( 3.17 ) for a massless scalar in AdS2- The 
in-going solution to eq. (|3.17] ) is 

c|(p) = e-f = l-!5 + of5y n = uVl-e^, (4.1) 



where in the second equality we have expanded in O/p ^ 1. The leading terms are a 
constant and p~^ , precisely as we expect for a massless scalar in AdS2i dual to a dimension- 
one operator. Following the standard AdS/CFT recipe, we identify the coefficient of the 
leading, non-normalizable term, the constant, as the source for the dual operator, and 
the coefficient of the sub- leading normalizable term, p"^, as the one-point function. The 
dual two-point function is thus the ratio of these coefficients, — ifi, which is the expected 
result [16]. The full inner-region solution for <&(/?), including corrections, is 

$(p) = pe-f (l + 0(p^)). (4.2) 
At low frequency we can approximate the solution as 

^{p)^p-iQ. (4.3) 



In the outer region we obtain the leading-order small-frequency solution from eq. (3.13a) 
by setting w = 0. The resulting equation will have two linearly independent solutions which 
we will denote $o(/o) and $i(/o). These will depend on the parameters (k,e,m) but obvi- 
ously not on 0,. We then construct the general solution by taking a linear combination of 
^o{p) and ^i{p), with coefficients that depend not only on (k,e,m), but also on Q. To 
match to the inner region, we demand that $o(p) approach a constant, say one, as p — >• 0, 
and that $i(p) approach p as p — >• 0. That fixes the 0-dependence of the ratio of the 
coefficients. Explicitly, we find at low frequency 

^p) = ^i{p)-inMp)- (4.4) 

Near the AdS^ boundary we will find two solutions, one normalizable and one not. We 
denote these ^v{p) and $5(/o), respectively, where "V" and "S" stand for "vacuum expec- 
tation value" and "source." These will be linearly related to <l*i(p) and ^o{p)- 

$i(p) = ai$s(p) + bi<^v{p), (4.5a) 
$o(p) = ao<l>sip) + bo<^>vip), (4.5b) 

where we have suppressed the (A;, e, m)-dependence of the coefficients ai, 6i, oq, bo. Near 
the AdS^ boundary we thus have 

^>(p) = (ai - inao)<^s{p) + {bi - inbo)<i>v (p) ■ (4.6) 
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The coefficient of the second term is the one-point function and the coefficient of the first 
term is the source, hence the retarded Green's function, G{Q), is simply the ratio 



G{n) = h.Z^, (4.7) 
ai — ulaQ 

Notice that we have suppressed G(J7)'s dependence on (k,e,m). Eq. ( p7| ) is the leading 
low-frequency form of the retarded Green's function for all operators dual to fiuctuations 
obeying eq. ( p. 13a ). Clearly the Green's function of the (O-l-l)-dimensional CFT, —iQ, 



determines the low-frequency behavior of the full Green's function. 

For <I>^(p), the analysis is nearly identical. The main difference is that in the inner 
region we do not extract a factor of p from <I>^(p), so eq. (4.1) is the inner region solution 
rather than eq. ( |4.3D , or in other words in the inner region, p and a constant are replaced 
by a constant and p~^, respectively. We denote the outer region solutions that match onto 
these as ^^(p) and ^^^{p). The rest of the procedure is unchanged. In the outer region 
$^(p) and ^^i{p) are linearly related to the normalizable and non-normalizable solutions 
as (in hopefully obvious notation) 

^^ip) = a^^^ip) + b^^${p), (4.8a) 
«>±i(p) = a%<^^ip)+b%^$ip). (4.8b) 

and the retarded Green's function at low frequencies is 

= I . (4.9) 

The coefficients ai, oq, bi, bo and Oq , a^^, 6q , b^^ are necessarily real, as is clear from 
the zero- frequency form of eqs. ( p. 13a ) and ( 3.13b| ). To obtain them requires solving for 



$(/o) and $(p)^ for all p, which we have only been able to do numerically. Crucially, 
however, to solve for the coefficients we can set Q = 0, in which case p = is a regular 
singular point of eqs. ( |3.13a| ) and ( p.l3b| ), for which we need no special methods to solve 
the equations. We may thus straightforwardly compute these coefficients for various values 
of (k,e,m) and, using eqs. (4.7) and ( |4.9D , we may look for instabilities. A QNM with 



negative imaginary part would appear when ^ 0" or ^ 0". 

We have studied numerically the m = 0, 1 modes of ^{p), the m = 1 mode of $"*"(/?), 
and the m = 1,2,3 modes of $~(p) for values of e G (0,0.9) and k G (0,50). For ^{p) 
and $"*"(/)), we detected no instabilities. For <&^(/o) we found that — )• D{e)k^ > 

as fc — >• 0, for some function D{e) that we can compute numerically (see figure |^). Such 
behavior does not indicate an unstable mode, but rather a mode with a dispersion relation 
Q = iD{e)W' . This is our R-spin diffusion mode. The operator dual to $~(p) with m = 1 is 
a dimension-two Lorentz scalar in the vector representation of SU (2)/j [29]. To be explicit, 
let = (g, q^)^ denote the SU{2)r doublet of squarks. '^~{p) with m = 1 is then dual to 

= Q^Vi^g^ (4.10) 
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where a^^ are the SU{2)r Pauh matrices. 

We mostly consider modes with m < 2. We can make a heuristic argument for why 
modes with higher m are unhkely to exhibit instabihties. Recall from section § that the 
value of m determines the dimension of the dual operator. Larger m means larger operator 
dimension. From a field theory point of view, we expect instabilities to appear first in 
the most relevant operators, i.e. the lowest-dimension operators, as we cool the system 
through a phase transition, for example. If we see no instabilities for small m, then we 
have good reason to expect that no instabilities will appear for larger m. 

4.2 The Zig-Zag Method 

We now come to the problem of computing the QNM's of our system numerically. The 
main obstacle is the presence of irregular singular points in eqs. ( ^.ISa ) and ( 3.13b ). Li 



particular, /O = is a simple pole of the coefficients pi{p) and P2{p) but a fourth-order pole 
of qi{p,U},k,m) and q2{p,u},k,m) 

Let us illustrate the problem using the near- horizon equation for ^{p), eq. ( 3.16| ). 



That equation has two solutions, /oe^'^/^, where the minus sign gives the in-going solution. 

If we introduce a cutoff at some small but finite p, impose the in-going wave boundary 

condition, and numerically integrate the equation of motion, the results are extremely 

sensitive to changes in the cutoff, and are thus unphysical. Moreover, since the QNM 

usually contains a non-vanishing imaginary part, the in-going solution is exponentially 

growing as p — )• while the out-going solution is exponentially suppressed. Since our 
_ — — 

series solution pe p {l + 0{p )) is only accurate sufficiently close to the horizon, numerical 
error will always effectively "include" the exponentially suppressed solution, leading to an 
inaccurate result. 

A few techniques have been developed to deal with such problems [16,50,51,60,79]. 
One is to use matched asymptotic expansions at low frequency, which we employed above. 
Another technique, valid for any frequency, is the matrix method, or Leaver's method 
[50,51]. Consider $(/>) as an example. With Leaver's method we begin by extracting 
from ^{p) the leading behavior near each singular point, so we set $(p) = pe-'^/P^{p). 
Inserting that into eq. ( ^.13a| ), we obtain an equation for ^{p). To solve that equation, we 
choose a non-singular point of the equation, p, such that when we expand ^{p) about that 
point the radius of convergence includes the two singular points at the horizon and at the 
boundary, 

n 

m = J2'^p(p-pr^ (4-11) 

p=0 

with some coefficients Cp that depend on u, k, e, and m. We truncate the series at some 
order n and then insert eq. (4.11) into the equation of motion. We obtain n + 1 equations 



for the n + 1 unknown coefficients, which we can assemble into a matrix equation, 

n 

X]-^gpCj, = 0, (4.12) 

p=0 
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where Aiqp also depends on u, k, e, and m. We then find the QNM's by finding where 
detAi = 0. Upon increasing n, the results should converge to the physical values. 

Leaver's method works well for fields in extremal AdS-RN [20,51,80]. For example, for 
a charged scalar field [51], the coefficients Cp obey a nine-term recurrence relation, so that 
Aiqp only has nonzero entries along a diagonal band of width nine. A numerical calculation 
of det is then reasonably efficient, even for large values of n. 

Leaver's method is more difficult to use in our case. In general, we found that Ai 
is upper triangular, indicating that the Cp obey a more complicated recurrence relation 
compared to an AdS-RN charged scalar. The calculation of M becomes more and more 
difficult as we increase n, since we have to compute the determinant of a large upper 
triangular matrix rather than of a reasonably sparse matrix with entries only along a 
diagonal band. 

To make progress, we developed our own method, which works as follows. First we pick 
some value O. We want to test whether that value is a QNM, so we denote it ^^qnm- Next 
we analytically continue p to complex values, such that argp = argcDq^^j^. The benefit of 
doing so is that in the exponent of pe^^^^^, the ratio fi/p is now purely real, so the in-going 
and out-going solutions are purely oscillatory and thus easy to disentangle. Now we just 
need to specify a contour for the p integration that takes us from p = to the point at 
infinity on the real p axis, starting along the line argiDq^j^ for small values of 

Complexifying p comes at a cost, however, since branch cuts appear in the complex p 
plane. These arise from the various factors of y{p) in the coefficients P2{p), Qi{p,ijj,k,m), 
and q2{p,uj,k,m). To see these, recall from eq. ( p. 12 ) that y{p) is an incomplete Beta 



function, which can be expressed in terms of a hyper geometric function. In our case, for 
real p we have 

B(J^^,^,)-,^^.F.(lll■.^). (4.3) 



p6 + l'6'3y (l + p6)i/6^ ^\^6'3'6'i + p6 

The hypergeometric function has a branch cut for values of its argument from one to 
infinity. Upon complexifying p, in our case these points correspond to p = c« and the 
points -|- 1 = 0, respectively. We thus find six branch cuts in the complex p plane, as 
illustrated in figure ||. 

The integration of the equations of motion in the complex p plane will obviously 
encounter problems if p crosses these branch cuts. We thus employ the "zig-zag" integration 
contour illustrated in figure |l[ For small \p\ we move along the line argWQ^*,^^. For some |p| 
large enough that the contamination from the out-going solution is negligible, but smaller 
than the point where the branch cut begins, \p\ = 1, we turn toward the real axis, moving 
along an arc at the fixed radius \p\, until we reach arg p = 0. We then continue integrating 
along the real axis to p = oo, or actually to some large cutoff value. 

We must specify matching conditions at each of the turning points of the contour. Let 
us denote the two turning points as pi and p2, and the value of $(/o) on the three pieces 
of the contour as ^i{p), ^ii{p), and ^ni{p) (see figure ||). To guarantee that solutions 
will be analytic in the complex p plane (away from the branch cuts), we require that the 
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Figure 1: (Color online) The complex p plane, with the branch cuts in the coefBcients of eqs. ( 3.13a ) 
and (3.13b) indicated by dashed blue lines and our contour of integration indicated by the solid 
orange curve, with segments I, II and III labeled, as well as the turning points pi and p2. The 
dotted black line ending in the dot labeled Wqnm indicates argWqg'^M- 



function and its first derivative agree at the turning points: 

^i{pi) = ^iiiPi), $'/(/0i) = ^'iiiPi). (4.14a) 
^ii{p2) = '^iii{p2), ^//(/02) = ^m(/02), (4.14b) 

where the derivatives are taken along the corresponding piece of the contour. 

To illustrate our procedure from start to finish, consider for example ^{p)- We first 
expand ^{p) near the horizon. 

Hp) = 6-'''/^ + 

where we calculate the constant coefficients recursively from up to ap^^^, with Pmax 
typically of order twelve. For the numerical integration, we typically imposed a small- 
p cutoff of p < 0.1, although the results of the zig-zag method were extremely stable 
against changes in the cutoff. We then integrate the equation numerically along the zig- 
zag contour up to some large-p cutoff, near the AdS^ boundary. We then fit the solution 
to the leading normalizable and non-normalizable solutions in the large-p region (explicitly 
restoring $(p)'s dependence on Q and k for the moment), 

Hp,c^,k) - A{Ld,k)^sip) + B{cd,k)^v{p). (4.15) 
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For example, when m = O/^ 

^sip) = l+^^^^^^^ + Oip-'\ ^v{p) = r' + Oip-'). (4.16) 

The coefficient A{uj, k) acts as the source for the dual operator, while B{uj, k) is the one- 
point function. The retarded Green's function, in the regime of linear response, is thus 
(proportional to) the ratio B(LD,k)/A{ijj,k). To locate the QNM's we looked for zeroes 
of \A{u!,k)/B{u!,k)\. We performed the search for the zeroes numerically using our own 
two-dimensional (over the complex u) plane) minimization algorithm.^'' Actually, in a 
numerical procedure \A{u}, k) / B{u),k)\ will never be exactly zero, so we searched for values 
of \A{Cj^k) / B{u,k)\ below some cutoff typically between 10~^ and 10~^. Our procedure 
for $^(p) was very similar to that of ^{p)- 

As mentioned above, in all cases the results of zig-zag method were extremely stable 
against changes in the small-/) cutoff. As a more rigorous test, we compared the results 
of the zig-zag method to the low-frequency forms of the two-point functions derived in 
section [4.1| . Since we only computed QNM's, and not the full two-point functions, we 
compared only to G~(0) in eq. (|4.9D because that has a QNM at low-frequency, the R-spin 
diffusion mode with dispersion relation = iD{e)k'^ . The results of the zig-zag method, in 
the low- frequency limit, produced precisely this dispersion relation for the R-spin diffusion 
mode. 

4.3 Numerical Results 

For ^{p) and m = 0, the first three QNM's appear in figure for e up to order one and k up 
to order ten. For larger k we were unable to identify zeroes of |j4((D, k)/B{co, k)\ numerically 
because the value of \A{u}, k)/B(Q, k)\ in the immediate neighborhood of the zero is larger 
by orders of magnitude, i.e. numerically we find a large "background" surrounding the 
zero, so that identifying the exact location of the zero is difficult. Notice that the QNM's 
are mirror-symmetric about the imaginary a) axis. 

None of the modes that we see ever cross into the lower half of the complex w plane, 
that is, none become unstable. One unusual feature we do see from figure ^ is that the 
imaginary part of the QNM decreases as the momentum k increases. (Contrast that with 



the zero sound dispersion relation in eq. (|2.1lD .) Figure y shows the real part of the first 



QNM as function of A;, for e G (0,0.9). At large k we see that the modes are approaching 
the lightcone. Numerically we find that by A; ~ 10 their slopes are within about 10% of 
the speed of light. We expect them to reach the speed of light in the limit A; — ?■ oo. 

We also studied the QNM's of ^{p) with m = 1. The results were qualitatively similar 
to those for m = 0, so we will not present them. In particular, none of the modes became 
unstable for the values of e and k we explored. 

^®For general m, at leading order $s(p) ~ and $v(p) ~ p^™"^. The leading behavior of ^v{p) 
agrees with that obtained for the scalar and type II and III gauge fluctuations in ref. [29], as expected. For 
$±(p), at leading order $+{p) = p'"+i and <I>+(p) = p-'"-^ while $s(p) = p""-^ and $^(p) = p-"'-^ 
The leading behaviors of $y (p) agree with those of the type I fluctuations in ref. [29] . 

'^'^A Mathematica notebook implementing the minimization algorithm is available from J. Shock upon 
request. 
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Figure 2: (Color online) The complex ui plane, where the curves show the motion of the first three 
QNM's of for m = 0, labeled by n = 1, 2, 3. The solid lines describe the motion as a function 
of momentum k G (0, 10) while the ten lines for each n (nine for n = 2) are for e — M/fi = to 
0.9 (0.8 for n = 2) in steps of 0.1. The dashed lines indicate the behavior for fc = as we increase 
e. None of the modes cross into the lower half-plane for our ranges of s and k, hence none become 
unstable. 




Figure 3: (Color online) The real part of the first quasi-normal frequency of ^{p) for m — (the 
point n = 1 in figure ^) as a function of k, for e = Af//i G (0, 0.9) in steps of 0.1 from bottom to top. 
The solid blue line intersecting the origin is the lightcone, Re(d)) — k. We see that as k increases 
the modes approach the lightcone. Numerically we find that at A: = 10 the asymptotic velocities 
are typically within about 10% of the speed of light. 



Our results for the QNM's of ^"^(p) with m = 1 are very similar to those for ^{p) with 
m = 0, although numerically the QNM's were more difficult to locate so we only present 
results for the first, in figure Clearly the first QNM never crosses into the lower half of 
the O plane and hence never becomes unstable. On general grounds, we thus do not expect 
the higher m = 1 QNM's, or the QNM's of ^~^{p) with higher m, to become unstable. 

Our results for the first three QNM's of ^~{p) with m = 1 appear in figure ^. The 
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Figure 4: (Color online) The complex ui plane, where the curves show the motion of the first 
QNM of <^+ip) with TO = 1 for e = (0,0.2,0.4,0.6,0.8) and k £ (0,10). Our labeling conventions 
are the same as in figure |[ The first QNM never crosses the real axis for the values of e and k 
that we consider, and hence never becomes unstable. The higher QNM's are more difficult to locate 
numerically. Given that the lowest mode never becomes unstable, and the higher modes are farther 
from the real ui axis, we do not expect them to become unstable cither, for our ranges of e and k. 
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Figure 5: (Color online) The complex Hi plane, where the curves show the motion of the first three 
QNM's of ^~{p) with m = 1 as we vary e and k. The lines have the same interpretation as those 
in figure ||, although in this case for the n — I and n = 3 modes we only went up to e = 0.8 and 
for the n = 2 mode we only went up to e = 0.7. We considered k e (0, 10). The behavior of the 
QNM's is similar to that of the $(/?) modes shown in figure ||. 

results are qualitatively similar to those for the <5(p) QNM's. In particular, none of the 
modes become unstable in the range of momenta we studied. The limiting velocities are 
qualitatively similar to those of the $(/)), m = QNM's (see figure ^). 

A special QNM appears for $~(p), qualitatively different from the others in that it 
sits precisely on the imaginary axis at small k. The motion of this mode as we increase k 
is illustrated in figure ^. For sufficiently small k, we have found via numerical fitting that 
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Figure 6: (Color online) The complex uj plane, where each line indicates the motion, as a function 
of momentum fc e (0,9), of the $~(p), m = 1 QNM that starts on the imaginary axis at small fc, 
our "R-spin diffusion" mode. The different lines are for e = A///i = (0,0.1,0.3,0.5,0.7,0.8,0.9). 



the dispersion relation is approximately 



w = i^M=P + /(e)P + C'(fc^), (4.17) 



e 



for functions D{e) and /(e) that we can determine numerically. The leading behavior is 
similar to that of a diffusion mode for a conserved current, with diffusion constant f ^^l. 

Vl— £ 

(in terms of our dimensionless frequency and momentum a) and k\ Notice, however, that 
<!>"(/>) with m = 1 is dual to the scalar operator in eq. ( 4.10| ) [29]. A similar mode, 
with the same leading-order dispersion relation, also appears among the QNM's of ^~{p) 
with m = 2,3. We plot the value of ^^^^^ as a function of e G (0, 1) for the three cases 
m = 1, 2, 3 in figure |^ 

We discuss the physical interpretation of this "diffusion mode" in the next section. 



5. Discussion and Conclusions 

We have studied the stability of holographic quantum critical matter in a top-down model, 
a single probe D7-brane in AdS^ x with nontrivial worldvolume gauge field, dual to 
= 4 SYM theory coupled to a single flavor of A/" = 2 hypermultiplet, in a state with 
zero temperature and finite baryon density. We computed the spectrum of fiuctuations 
of D7-brane fields, dual to the spectrum of mesons. For the range of momenta that we 
studied, k G (0,50), we detected no instabilities. 

Beyond that, we had three main results. The first was to show that all bosonic world- 
volume fluctuations effectively "see" a near-horizon AdS2 , due to the nontrivial background 
worldvolume flelds. The dual field theory statement is that in these finite-density states, 
some (0-|-l) dimensional CFT emerges at low energies. The second main result was our nu- 
merical procedure, the zig-zag method, which may be useful for any calculation of QNM's 
for a bulk differential operator with singular points. Our third main result was the detec- 
tion of a purely imaginary mode, similar to the diffusive mode of a conserved current, in 
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Figure 7: (Color online) The "diffusion constant" defined in eq. (i.l7), coming from the 

m = 1, 2, 3 QNM's that lie on the imaginary axis at small momentum fc, versus e = Mj p.. 



the two-point function of the operator in eq. ( [4.10 ), a dimension-two scalar operator that 



is a vector of SU{2)r. If we interpret SU{2)r as spin, in analogy with electronic systems 
at low energy, then the mode we found has a natural interpretation as the diffusion of 
the "R-spin." For higher values of m, the dual operators form an (m -|- 2)-dimensional 
representation of 5f7(2)/j. For all of m = 1,2,3, we saw a diffusion mode (see figure |^, 
and we generically expect such a mode for all m. 

An important task for the future is to extend our analysis to Nf > 1. The flavor 
symmetry is then U{1)b x SU{Nf), so we have the option of introducing a chemical po- 
tential for each Cartan generator of SU{Nf), which in QCD language would be isospin 
chemical potentials. Here we expect that at zero temperature the ground state should in 
fact be a p-wave superfluid, representing a condensate of vector mesons [74-77]. Some 
natural questions then are whether any scale invariance emerges at low energies, i.e. a 
(O-l-l)-dimensional CFT, and whether an R-spin diffusion mode appears.^® 

Another important task for the future is to extend our analysis to other types of probe 
D-branes. An especially attractive example is a probe D5-brane extended along AdS^ x 5^ 
inside AdS^ x [81-83], because an exact solution for the worldvolume fields, similar to 
the solution in eq. ( |2.8D , is known in that case [28] . Here the dual flavor fields are restricted 
to a (2-|-l)-dimensional defect. The spectrum includes a zero sound mode [30], and at very 
low energies we expect to find an emergent (O-l-l)-dimensional CFT [33]. The spectrum of 
fluctuations will obviously be different from that of the D7-brane, and instabilities could 
appear. In particular, an instability involving a violation of the AdS2 BF bound may occur, 
as happens for the D5-brane with an external U{1)b magnetic field [33,84]. 

Our results add to the list of properties characterizing whatever microscopic state the 
solution in eq. ( ^.8)) represents. That state has a low-temperature heat capacity cy oc 



^*For Nf — 2 and at finite temperature, an isospin chemical potential, if represented in the bulk by a 
constant gauge field on the D7-branes, produces a instability towards squark condensation [68]. Whether 
the same instability appears for a radially-varying gauge field remains an open question. 
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T^. The low -energy spectrum includes not only the zero sound mode but also the R- 
spin diffusion modes. At the lowest energies, some (0+l)-dimensional CFT emerges which 
controls the low- frequency form of all correlation functions. Despite having a nonzero 
entropy at zero temperature, the state is stable against small perturbations. 

What effective theory replaces hydrodynamics at low temperatures? That theory must 
include not only the zero sound mode but also the R-spin diffusion modes. Can a Bose- 
Einstein condensate of strongly-coupled squarks, or some other more exotic microscopic 
physics, give rise to the properties listed above? In particular, what is the microscopic 
origin of the (0+l)-dimensional CFT? We plan to study these and related questions in the 
future. 
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Appendix: the Zero Sound Mode 

In this appendix we re-derive the results of refs. [30,31,53] for the dispersion relation of the 
zero sound mode, which among other things provides a consistency check for our fluctuation 
equations of motion, eqs. ( 3.6a ) to ( 3.6| ). 



In the field theory, the zero sound mode appears in the U{1)b current two-point func- 
tions. These operators are singlets of the 50(4) global symmetry. They correspond to 
fluctuations 5 At and (5^3 that are zero modes on the 5^, that is, fluctuations with m = 0. 
When the flavor mass is nonzero, so that the D7-brane embedding is non-trivial, these 
fluctuations couple to the m = mode of the fluctuation x as well [31]. This suggests that 
the ansatz in our eq. ( ^.8a ) should include the zero sound mode, however that is not the 



case. When m = in that ansatz, eq. (3.6f) becomes trivial since the derivatives acting in 
directions vanish. Eq. (3.61) is essentially a constraint equation that comes from setting 
5Ai = 0. When a constraint is trivially satisfied, we expect to find an additional physical 
degree of freedom. In other words, an m = fluctuation distinct from that of eq. ( 3.8a| ) 
should appear. 

We thus consider an ansatz in which only the fluctuations X) ^At, 5A^ are nonzero, 
all are zero- modes on the S^, and all are plane waves of the form g*'^"*'^^^ in Minkowski 
directions. In this appendix, for all functions we will suppress all dependence on the 
coordinates /j, t, and X3 unless otherwise noted. 
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We will work in terms of two functions, /i and /2, defined as 



/i = 5 At - ex, f2 = cor^'^'y'x + '^rf'lrf'^At - UA^. 
In terms of these, the constraint equation, eq. ( |3.6c ), becomes 

5p/2 = a;5,(7?*Vr?--)(/i). 
The relevant equations of motion, eqs. ( |3.6b| ) and ( 3.6d| ), may be written as 



(1) 



(2) 



flb XX pt I ' 



xx2 pt / 



= 0, 
(3) 

da{V^r,fr,-''r,P'y'dtx) + daiV^VsV''''{v''/vndbSAt) + u:^rf'i^^- h = 0- (4) 

^ [(3) + e(4)], and making use 



Taking the linear combination of the above equations 
of the identities 



£ 



-11 + ^ 



we obtain /I's equation of motion, 



da 



0. 



(5) 



(6) 



Crucially, notice that our fields fi and /2 are not gauge- invariant, therefore we must also 
consider pure-gauge solutions [85]. The advantage of our approach is that we only have 
one dynamical equation, eq. (^), in contrast to the approach of ref. [31]. 

Notice that when the flavor mass M = 0, so that e = 0, x disappears from eq. (|^). As 
explained in section ^, that occurs because the chiral U{l)ji symmetry is restored when 
M = 0, and x is charged under J7(1)_r but 6At and SA3, are not. In the M = limit, 
we can identify the above fluctuations of the gauge field, along with 6A1 and 6A2, as the 
finite-density generalization of the type II fluctuations of ref. [29], with m = 0. 

To find the zero sound dispersion relation we follow refs. [30,31]: we solve the equations 
in two different limits and then match the solutions in a regime where the limits overlap. 
More precisely, we first take a near-horizon limit of the equation, — >■ 0, and then expand 
the solution for low frequency and momentum, meaning w/p <C 1 and fc/p <C 1 with uj/k 
fixed. We next take the low-frequency and low-momentum limit of the equation, solve that 
for all p, and then expand the solution near the horizon to perform the matching. 



In the near-horizon limit, the in-going solutions are, with = wVl — e^, 



/i = ci e-f p'{l + 0{p)), 



f2 = C2e-- p\l + 0{p)), 



(7) 



where the equations of motion dictate that the normalization constants ci and C2 are related 
as 

6x 



C2 



(l-e2)i/2e2^2 



Cl, 



(8) 
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where N is defined above eq. ( 3.11 ). We expect to find such a relation, since we have only 
one dynamical equation, and hence only one integration constant. 
In the low- frequency and low- momentum limit, eq. (^) becomes 

{e^ - l)dM' hn^/^dph + dpi^^'ii^'^dph) = 0, (9) 

for which the solution is, in terms of the dimensionless radial coordinate p defined in 
eq. CT . 

where C is an undetermined constant. Notice that eq. (^) depends only on so we can 
add any constant to /i to obtain a new solution. We have fixed the constant by demanding 
that the near- horizon limit of the solution match onto the low- frequency limit of the solution 
m eq. (0), that is, we demand that the two solutions agree in the regime where the two 
limits overlap, as advertised. The integral in eq. ([l^) describes an incomplete Beta function, 
although we will not need that representation here. Integrating the constraint equation 
gives /2, 

fp 6 _^ _^ 

h{p)=^ ^Tttt -^fiip)dp. (11) 

Jo P \^ ~ ^ ) 

fi and /2 approach constant values at the AdS^ boundary, 

r(l)r(l) 2w r(^)r(|) ^ ^ 

lim/i = C^^^, lim/2 = C- 



As emphasized earlier, we must also consider pure gauge solutions since our fields 
are gauge-dependent, and change under gauge transformations 5Af, — t- 6Ah + dbO. To 
preserve our gauge choice 6Ap = 0, a must be of the plane- wave form e*'^*~*'^^3^ with no p 
dependence. The pure-gauge solution is thus 

X = 0, 5At = ioje'^'~'^='\ = -i/te^'^*-^^^^ (13) 

Recalling the definitions of /i and /2 in eq. ([l|), for the pure-gauge solution we find near 
the AdS^ boundary (where rf^ jrf'^ — )• —1 in eq. (|ip) 

lim /i = iw, lim = -i(w^ - k^). (14) 

The leading behavior of f\ and /2 near the boundary will be a linear combination of 
eqs. (p^) and (|l^). These leading constants act as sources for the dual operators, which 
roughly speaking will vanish when cj is a QNM. More precisely, when a; is a QNM, the 
following condition holds: 
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which gives us the leading term in the dispersion relation, 



1 



3 - £2 - 3^2 _ JVf 2 ^ 



(16) 



which agrees with the results for the speed of zero sound squared in ref. [31] when M 7^ 
and in ref. [30] (and eq. (pl|)) when M = 0. 

We should mention that, like all of the other bosonic worldvolume fluctuations, the 
fluctuations giving rise to zero sound "see" an effective AdS2 near-horizon region [21]. That 
is clear from eq. (0), since the e~*^/^ factor is a solution of eq. ( |3.17D , indicating that after 
removing a power of the radial coordinate, in /I's case and in /2's case, in the near- 
horizon region the fluctuations obey an equation of motion identical to a massless scalar 
in AdS2- 
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